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Abstract 

In this paper we study a reflected Markov-modulated Brownian motion with a two sided 
reflection in which the drift, diffusion coefficient and the two boundaries are (jointly) 
modulated by a finite state space irreducible continuous time Markov chain. The goal 
is to compute the stationary distribution of this Markov process, which in addition to 
the complication of having a stochastic boundary can also include jumps at state change 
epochs of the underlying Markov chain because of the boundary changes. We give the 
general theory and then specialize to the case where the underlying Markov chain has 
two states. Moreover, motivated by an application of optimal dividend strategies, we 
consider the case where the lower barrier is zero and the upper barrier is subject to 
control. In this case we generalized earlier results from the case of a reflected Brownian 
motion to the Markov modulated case. 

Keywords: Markov modulation, Brownian motion, dividend payout, two sided 
reflection 



1. Introduction 

A double sided reflected process, say at zero from below and some positive level b, is 
a reasonable model for a storage process where the stored quantity has to be nonnegative 
and the buffer size is limited. When borrowing or backlogging is allowed, then the lower 
barrier could also be negative. There is a huge literature on such processes, in particular 
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when the driving process (before reflection) is Brownian motion. Less attention is given 
to the case where the boundaries are themselves stochastic processes. For most papers 
on this topic the focus was on showing the existence and uniqueness of solution of the 
related Skorohod problem. A recent study which refers to many of the earlier results in 
this particular direction is [10] where the focus is on multidimensional models. For the 
one dimensional double sided reflection (non-modulated case) , we mention the important 
results reported in [9] and references therein. 

Very little work is done related to the computation of the stationary distribution of 
such processes when more explicit stochastic structure is assumed, especially when the 
boundaries are not smooth. One example of such a study is given in [8] where the driving 
process is Levy and there is only one lower boundary which increases linearly and then 
drops back to zero at arrival epochs of a Poisson process. 

We are not aware of any results for the case where the boundary together with the 
driving process are jointly modulated by some other process. This is motivated by 
situations in which the buffer size and the allowed backlog are allowed to change from 
time to time as a response to changes in the driving process which are caused by changes 
in an underlying environment. 

In this paper we model the environment as a finite state space irreducible continuous 
time Markov chain. When in a given state, our process behaves like a two sided reflected 
Brownian motion with drift and diffusion coefficient as well as lower and upper boundaries 
which are allowed to depend on this state. The main goal is to give a computational 
scheme for computing the joint stationary distribution of the buffer content and the state 
of the underlying environment. 

The paper is organized as follows. In Section 2 we present the general model and 
provide some preliminary results. Section 3 is about the stationary joint distribution of 
the buffer content and the underlying environment. Section 4 specialized the results to 
various cases where the lower barrier is zero (no backlog) and the underlying environment 
changes between two states. Under some conditions, for this case we also show how to 
compute the distribution of some regenerative epoch associated with this process. Finally, 
in Section 5 we generalize results of [6] , who considered the upper barrier as a cutoff point 
above which a company must pay dividends that are modeled by the regulating process 
at this upper barrier. 

2. Model 

Let W — {W(t)\ t > 0} and J = {J(t)| t > 0} be two independent processes 
where W is a Wiener process (a standard Brownian motion) and J is an irreducible 
and homogeneous continuous time Markov chain with state space E = {1, . . . , N}. We 
assume that J has right continuous sample paths and we denote by Q = (g,j) its rate 
transition matrix, by 7? = (7Tj) its stationary distribution and define P = diag[7?]. For 
each i G E we let a(i) < b(i) be two finite real numbers which define the upper and lower 
barriers when in state i. 

It is standard to show that there is a unique process (Z, L, U) satisfying 

Z(t) = Z(0) + f a{J{s))dW{s)+ f n{J{s))ds + L(t)-U(t) (1) 
Jo Jo 
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where a(J(t)) < Z(t) < b(J(t)) for each t > 0, L and U arc nondecreasing right contin- 
uous processes with L(0— ) = U(0—) = 0, 

J (z(s)-a(J(s)j)dL(s)=0 and J (b{J{s)) - Z(s)) dU(s) = . (2) 

Z is the two-sided (Skorohod) reflection of the modulated process 

X(t) = Z(0)+ [ a(J{s))dW(s)+ [ n(J(s))ds (3) 
Jo Jo 

at the modulated barriers [a(i), b(i)], i G E. We denote by k = the asymptotic 

drift of the process X{t). 

Although Z is not Markovian, (Z, J) is. Let us identify its generator. 

Let f(w,i) be a bounded twice continuously differentiable function in w satisfying 
f'(a(i),i) — f'(b(i),i) — 0. We note that clearly there always is a twice continuously 
differentiable h : M 2 -> M such that h(w,i) = f(w,i). The generalized Ito formula 
for semimartingales (e.g. Theorem 33 on p. 8 of [11]) now implies after some obvious 
manipulations that 

f(Z(t),J(t)) = /(Z(0),J(0))+ /' a(J(s))f'(Z(s),J(s))dW(s) (4) 

Jo 

+ J* ^^-f"(Z(s), J(s)) + M (J( S ))/'(Z( S ), J(s)) ds 
+ f f'(Z(s),J(s))d(L c (s)-U c (s))+ A/(Z( S ),J( S )) 

J ° 0<s<t 

where for cadlag functions of bounded variation on compact sets g we denote g(s—) = 
lim uts ff(u), Afii(s) = g(s)-g(s-) and 5 c (s) = g(s) - Eo<«< ;s A ff( u )- Now > observe that 
f'(a(i), i) = and the first relation in (2) imply that 

* f(Z(s), J(s)) dL c (s) = [ f(a(J(s)), J(s)) dL c (s) = (5) 

Jo 

and similarly J * .f(Z(s), J(s)) d?7 c (s) = 0. 

Next, note that the only way that s can be a jump epoch of f(Z(-), </(•)) is if it is 
a jump epoch of J. If a(J(s)) < Z(s-) < b(J(s)) then clearly Z(s) = Z(s-). If cither 
Z(s-) < a(J(s)) or Z(s-) > b(J(s)), then in the first case necessarily Z(s) = a(J(s)) 
and in the second Z(s) = b(J(s)), therefore we can always write Z(s) = o(J(s))VZ(s-)A 
b(J(s)). This implies that 

Af(Z(s), J(s)) - f(Z(s-), J(s)) - f(Z(s-), J(s-)) , (6) 

where we used the notation f(w,i) — f(a(i) V w A b(i),i). Now for each i ^ j let Nij 
be independent Poisson process with rate g^, independent of W, such that if an arrival 
finds J in state i then it instructs J to jump to j and otherwise nothing happens. Now, 
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Nij(s) — qijS is a martingale with respect to the filtration generated by W and Nj for 
i ^ j and f(Z(s),i) are bounded processes. Hence, 

/ {f(Z(s-),j)-f(Z(s-),i))l {J{s _ )=i} d(N ij (s)-q ij s) (7) 

J[0.t] 

is a martingale and thus, by summing over i, j G E, with i 7^ j, and noting that = 
— J2j^i Qij > this implies that 

£ Af(Z(s),J(s))- f t J2<U(s-)jf(Z(s-),j)ds (8) 
o<s<t ^° jge 

is a martingale with respect to the filtration generated by W and the Poisson processes 
N^, but also with respect to the filtration generated by W and J, since it is adapted 
to that filtration and are non-anticipative with respect to it (see also [12], Lemma 
V.21.13). Since f'(w,i) are locally bounded, it follows that the following is a martingale 
with respect to the filtration generated by W and J: 

M(i) = Yl A/(Z(*),J( S ))- fY,<lJ(s)J{Z{s),j)ds 
o<s<t Jo jeE 

+ f a(J(s))f'(Z(s),J(s))dW(s) . 
Jo 

Rewriting equation (4) in terms of the above martingale we have 

f(Z(t),J(t))-f(Z(0),J(0)) = f Af(Z(s),J(s))ds + M(t) 

Jo 

where we denoted by Af(z,i) the following operator 

Af(z,i) = l<T 2 (i)f"(z,i) + n(i)f'(z,i) + Y,qijf(z,j), {z,i)&£. (9) 

j£E 

This gives the expression for the generator of (Z, J) restricted to the set of bounded 
functions twice continuously differentiable on the continuous component and satisfying 
f'(a(%),i) = f'(b(%),i)=0,i€E. 

We denote by £ — \J ieE [a(i), b(i)} x {i} the range of values assumed by the process 
(Z, J) and define the following subsets of the state space of J, 

E+ = {j e E; a(j) > or > 0} , and E~ = {j G E; a(j) > or < 0} . 

3. Stationary Distribution 

If « G E~ , a(i) < Z(0) < b(i) and J(0) = i, then the probability of hitting a(i) before 
J changes state is bounded below by the positive probability of this event starting from 
Z(0) = b(i) and J(0) = i. This observation together with a geometric retrial argument, 
recalling that J is irreducible, implies that [Z, J) is a regenerative process with finite 
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mean regeneration epochs. A similar argument can be made for i G E + . If E~ U E + = 
then the process is just a deterministic function of J and thus clearly positive recurrent. 
Thus, in any case a unique stationary distribution exists. 

We show that if a solution to (10) exists then it must be the (unique) stationary 
distribution. Later we will show how to construct it. 

Theorem 1. The stationary distribution of the process (Z, J) is the unique solution of 
the following system of differential equations 

\o\i) n'l{z) - n' t {z) + ]T qji n 3 (a(j) v z a b(j)) = o a(i) < z < b(i) (10) 

with boundary conditions ILj(a(i)) — and IL(6(i)) = ixi, i G E . 

Proof. The stationary distribution satisfies the following equation for any function / 
belonging to the domain of the generator A, 



Mi) 

Af(z,i)dll i (z) = , (11) 



Let / be any twice continuous differentiable function on £ with bounded support, then 
using integration by parts we get that equation (11) reduces to 



]T / ia 2 wn^)-M0n^)d/(^) + ]T / n^) df(z,j) = o . (12) 

Noting that 



/ 

J a 



Hi) . Mi) 

Il i (z)df(z,j)= Ui(a(i)VzAb(i))df(z,j) 

la(i) Ja(j) 

and interchanging the indexes i and j in the two sums in the last term of (12) we get 
that (10) implies (11). 



The way of solving the system (10) is to divide the interval [mini{a(i)}, maxi{6(i)}] 
to disjoint subintervals, to get a solution of the above system in each of them and then 
to appropriately glue together all these partial solutions. For this we set Iq = mirii{a(i)} 
and by lk+i = mini{a(i) > Ik} A mini{b(i) > Ik} and we define the closed intervals 
h = [h-uh], k = l,...,K. 

Fix one of these subintervals, say Ik, and let Ek = {i G E : a(i) < lk-i < h < b(i)]} 
be the set of states active over Ik- The restriction of the system (10) to the subinterval 
Ik then reads as follows 

^a 2 {i)U'f(z) - M(*)n^) + J2 1^ U j( z ) = c fcW fc-i <z<l k , i£Ek (13) 

jeE k 

where we set Ck(i) — J2j-b(j)<i k Qji w i- ^ ne equivalent matrix form of (13) is the following 

S k TL'l{z) - M k ti' k {z) + Q^n fe (z) - 4 (14) 
5 



where S k = diag[i £ E k : cr 2 («)/2], M k = diag[i e E k : Q k = G E k : q io )] and 

where Il k (z) and c k denote the restrictions of the vectors U(z) and c to the only states 
active over I k . [-] T is used to denote the transposition operator. 

If the set of active states over I k is a proper subset of E we have that Q is a strictly 
substochastic matrix and this has an inverse, therefore the system (14) admits the con- 
stant k k = [Q _1 ] T c as particular solution. In the case E k = E, Q reduces to the rate 
transition matrix of J that is stochastic and singular. For this case the constant c is zero, 
the system (14) is homogeneous and the particular solution k k = is the zero constant. 

Adding to the particular solution the homogeneous solution, we can always write the 
general solution of (14) in the following form, see [7], 

II fc (z) = T k e AkZ u k + k k . (15) 

where (T k ,A k ) is a Jordan pair of the matrix polynomial (14) and u k is the unknown 
vector that can be determined by the boundary conditions. For the specific case of second 
order matrix polynomial, a Jordan pair consists of a pair of matrix with the following 
properties 

S k T k A 2 k + M k T k A fe + Q fe r fe = o 

and the rank of col[Tfc, T k A k ] is maximum. 

To be able to glue together the solutions of all intervals I k , we will see later that it is 
necessary to be able to solve a system of equations, and this requires to prove that there 
exists the inverse of the matrix associated to the system. 

Let P k be the restriction of the matrix P to the states active on I k , the main 
ingredient follows from a slight generalization of the results in [4] where it is shown that 
the systems 

S k GL 2 TM k GL + P^QjP k G = (16) 

admit solutions (r^,A^) that are unique under the restriction that the matrices A ± 
are substochastic. In this case the matrices are transition matrices and in particular 
E ± T ± = I ± where the projection matrices E^ are defined as the submatrices of the 
identity matrix l k constructed by selecting only the rows corresponding to the states 
contained in the set E k D E ± . The substochastic nature of the matrices A^ is the special 
characteristic that we are going to use later to prove the non singularity in the final 
system of equation. 

The solutions (rj ,Aj) allow to immediately construct the Jordan Pair in (15) as 
follows 

r fc = [P fc r+ P k TH]; A fe = diag[A+,-A,7] . (17) 

This construction is always valid but in the case when A^ have in common the null 
eigenvalue. This happens only in the case E k = E and the asymptotic drift of the 
modulated process X{t) is zero, that is n = 0, see also Section 7 in [4]. We are going to 
exclude this case as it has no added difficulty if not the one of making more complex all 
formulas. 

Using the special selected Jordan pair in (17), the solution (15) over the interval I k 
can be rewritten in the following form 

U k (z) = P k T+e+ A i z u + + P k T^e- A * z uf + k k , (18) 
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and it is fully specified after assigning the unknown boundary values 



IIj(/ fe _i+) for any i G E k n E + and Ui(l k -) for any i £ E k Ci E~ . (19) 

To glue together the solution over the entire interval [IqJk] it is then necessary to 
solve for the unknown boundary conditions in (19) by using the constraints IF(a(i)) = 
and Tli(b(i)) = iri, for any i G E, together with the additional conditions on the continuity 
of the distribution functions 

Ili(aO')-) = Ui(a(j)+) a(i) < a(j) < b{i) (20a) 
IIi(60')-) = MKi)+) a(i)<b(j)<b(i) (20b) 

for any i G E + U E~ and j G E, and on the differentiability of the distribution functions 

WM.1)-) = n{(a(j)+) a(i)<a(j)<b(i) (21a) 
WM3)-) = Hj(6(j)+) a(i) < 6(j) < (21b) 
for any i G -E + fl E~ and j G E. 

3.1. Computing the stationary distribution 

In this section we show how to determine all the unknowns required to get the unique 
stationary distribution. We are going to consider the sequence of subintervals I k , with 
k = 1, . . . , K, and we represent the solution U k (z) in the fc-th interval according to the 
form (15). It is worth to say that when applying the model to specific cases it is possible 
to compute numerically the Jordan Pair in (15) or in its special form (17), like it is 
shown for example in [13] or in [1]. We used Wolfram Mathematica©, to get some of 
the analytical solutions shown in the examples of Sections 4 and 5, but for numerical 
computation any technical computing system, such as Matlab©, can be efficiently used 
for this purpose. 

In the same way as we defined the projection matrices we define the projection 
matrices for the states that belong to the intersection of the intervals I k -i and I k , that 
is the matrices D k , D\ and D k They are defined as the submatrices of the identity 
matrix 1^ constructed by selecting only the rows corresponding to the states contained 
respectively in the sets E k n E k -i, (E k D E k -i) n E+ and (E k n n (E+ n E~). 

For the states in E k that do not belong to the intersection E k _\ n E k we define the 
additional projection matrices D k and D k defined as the submatrices of the identity 
matrix 1^ constructed by selecting only the rows corresponding to the states contained 
respectively in the sets E k n Ek-i, (E k n -Efc-i) fl E + . 

Finally considering the intersection between the intervals I k and I k +\ we define the 
corresponding projection matrices U k , , U k , U k and U k whose definitions are easy 
to guess. In all definitions we have assumed that E = Ek+i = 0. 

Applying the continuity and differentiability constraints we get 

r>X r fe e ; *- iA * u k = Dt k k 

U k r fe e' fcAfc u k = U k (4 + 7T fc ) 

U k T k e l * A * u k -D k +i T k+1 e l * A ^ u k+1 =U k k k - D k+1 k k+1 

U k T k A k e lkAk u k -D k+1 T k+1 K k+1 e lkAh + 1 u k+1 =U k k k - D k+1 k k+1 

7 



for any k = 1,...,K, with 7?^ = col^, i G Ek], that defining 

A® = co\[Dt T k e l ^ Ak ,U^ T k e l * A \U k T k e l ^,U k T k K k e lkA % 
A u k = col[O k+1 ,O k+1 ,D k+1 T k+1 e hA ^\D k+1 T k+1 A k+1 e lkA ^], 
b k = col[D k k k , U k (k k + fffe), U k k k - Dk+ih+i, U k 0k], 
can be rewritten as 



A? u fe - u k+ i = b k ■ 



Defining the block matrix A whose block diagonal and upper block diagonal are 
made respectively of A k and A k with k = 1, . . . , K, all the unknowns may be solved by 
resolving the linear system Au — b, with u = co\{u k ,k = 1, . . . , K], and b = col[6fc, k = 
l,...,K]. 

The assumption that E + U E~ is not empty and that we can choose i E E + U E~ 
such that a(i) < b(i) implies that the matrix A has dimension at least one. In addition, 
considering the interval Ik such that a(i) = l k -\ if i e E + (corr. b(i) = l k if i S 
the corresponding fc-block in the matrix A contains a strictly positive submatrix 
£> k Tke lk ^ lAk (resp. t7 fc r^e' fcAfc ). Having that all other submatrices in the same rows 
are all zeros, we deduce that the matrix A is invcrtible if the rank of the fc-block column is 
maximum. It is then easy to realize by induction that A is invertible if for all k = 1 , . . . , K 
the corresponding fc-block column has maximum rank. 

Consider the k block, in order to prove that its rank is maximum we can look at 
a square submatrix with its dimension equals to the size of the block. Noticing that 
the matrices D k (rep. U k ) is a submatrix of D k (resp U k ) that corresponds to the 
complementary states in E + n Ek (resp. E~ C\E k ) not selected by D (resp. U ), after 
reordering the rows we get that the sought square submatrix is given by 



-iA+ 



E+PkTl 



\ E k P k T+e lkA i e- lkA * ) ' 

The matrix A k is invertible if its Schur complement is. It is given by 

g i fc -iA+ ^ _ e -h-iA.+ ^p^-^h-rA- e -l k A+ E -p kT + e -l k A~^ 

and it is invertible by the Levy-Desplanques theorem (see Lemma B.l in [14]) after 
noticing that all the matrices in the last term of the equation above are substochastic 
with at least one strictly substochastic. 



4. Fluid queues with modulated buffer 

One direct application of the model presented in the previous sections is the case of 
the fluid queue with Markov modulated buffer. Indeed if we assume that for all states 
the lower barrier is zero, i.e. a(i) = for any i e E, we can look at the process Z as the 
buffer content of a fluid queue whose net-input flows is the process X(t) defined in (3) 
and whose buffer is equal to b(i) when the environment J is in state i G E. 

In this simplified setting, it is easier to solve the system (10), as it is possible without 
loss of generality to order the buffer levels in increasing order, i.e. 6(1) < 6(2) . . . < 6(7V). 

8 



While the results of previous section are completely general and allow in any case 
to compute at least numerically the stationary distribution for any configuration of the 
barriers and the parameters of the Markov modulated Brownian motion, we limit the 
present discussion only to some few cases where it is possible to write the general solution 
in an easy form. 

It is worth mentioning that a closed form solution can be obtained for the general 
case when N = 2, but the expression is cumbersome and we prefer not to include it. It 
requires solutions of forth-order algebraic equations, sec Section 5.3 for an example of 
this type but applied to the dividend-payout problem. For this reason we further limit 
our focus to more specific cases and provide the solution for the case where the net-input 
flow is not modulated, i.e. the modulation applies only to the buffer levels, and for the 
case where only in one of the two states the diffusion coefficient is positive. 

4-1. Two sided Markov modulated reflection of a (ji, a) -Brownian motion 

In this section we analyze the case when E — 2 and the drift and diffusion components, 
\x and cr, do not depend on the environment. The reflecting levels are given by a(l) = 
a(2) = and < 6(1) < 6(2), and the system can be looked at as a fluid queue that 
for exponential periods of times, i.e. when J = 1, its buffer is reduced to a smaller 
value. A typical application for this could be a service station that for specified period 
of times receives help from another service station with which it shares the buffer size. 
When the second system turns on, the buffer of the first station is reduced and the 
overflow fluid becomes the starting content for the second station. In the following we 
derive the stationary distribution of the system, and in the next subsection we compute 
the discontinuity rate of the content process together with the size distribution of the 
discontinuities. 

The system (10) can be solved by considering the two intervals 1\ = [0, 6(1)] and 
I 2 = [6(1), 6(2)]. For the interval I\ we have that 

\o 2 n'{(z) -tin[{z) - qi2 nx(z) + q 21 n 2 (z) = o 

\a 2 n 2 '(z) -fiU' 2 (z) + q 12 IIi(z) - q 21 U 2 (z) = 
and for the interval I 2 we have 

Vn 2 '(z) -fin' 2 (z) -q 21 U 2 (z) = -g 12 7ri (23) 

In addition we have the boundary conditions 11,(0) = and IF(6(i)) = 7T; = qz-i,i/{qi2 + 
<?2i), i = 1)2, the continuity condition n 2 (6(l)— ) = n 2 (6(l)+) and the differentiability 
condition n 2 (6(l)-) = n 2 (6(l)+). 

Let Af = ©i ± Ai be matrix solutions of the equation \a 2 L 2 =F A* L + Q = ®- I R 
this case since = E^ = E we have that = = I, in addition Pi = diag[7ri, tt 2 ]. 
When /i < the matrices Af below are substochastic and are the same ones appearing 
in (18). When [i > they are not substochastic but the result still holds. When /i = 
it follows that the asymptotic drift k = and the root has multiplicity 2. As already 
mentioned before, for this case the solution looks slightly different from the one given 
below, and it can be computed using similar arguments. We have decided to omit its 
expression here. 
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(22) 



( 7T2 


-7T2 


MS 


7T2 \ 


V _7r i 


^1 , 




7T2 / 



Let A^ be the negative solutions of the equations |<7 2 z 2 ^/iz - (<?i2 + 92 1 ) =0. Define 
9i = (A+ + Af)/2 = -( A /// 2 + 2(g 12 + g 21 )(7 2 )/ < 7 2 and Ai = (A+ - A^)/2 = A = m/,7 2 , 
then we have that 

0! =6! 

and Ai = AI. Considering the continuity conditions at z — and z — 6(1) we can write 
the solution for z <G I\ as 

*g ) = e M^(i)) Pl sinh(0lZ ) s inh-(0 1 Kl))Pr 1 ( n^;^)) ) • (24) 

Let \ 2 be the negative solutions of the equations \a 2 z 2 =F \iz — q 2 \ = 0. Define 
©2 = (Xt + A 2~)/ 2 = -(vV + 2q 21 cr 2 )/a 2 and A 2 = (A+ - A^)/2 = A = [ija 2 , we can 
write the solution for z <G I 2 with boundary condition n 2 (6(2)) = ir 2 as 

n (-A * ^(z-b(i)) sinh(e 2 (6(2)-z)) 

H2(Z) = - 2 - e ( sinh(e 2 (6(2)-6(l))) (7r2 " n2(6(1))) ' (25) 

In the equation above we already used the continuity condition for z — 6(1) by imposing 
the boundary value at 6(1) to be equal to n 2 (6(l)). 

The solution is completely identified by solving for the value of n 2 (6(l)) that assures 
the differentiability of U 2 (z) at z — 6(1). 

We have that, for z G I\ 

U'(z)= AII(z) + e A ^- fc ( 1 »P 1 1 cosh(0 1 z) sinh- 1 (0 1 6(l))Pr 1 n(6(l)) , 
and for z e I 2 

DIM - - Afe - n, W) + e ^.->(.»e 2 -gMzl| y fe - n 2W i))) . 

To simplify the resulting expression define fc, = e 2 T Pi©i coth(0i 6(1))P^ 1 e?j, i — 1,2, 
where e?j is the i-th column vector of the canonical base in R 2 , and k 3 = 62 coth(62(6(2) — 
6(1))), then the value of n 2 (z) at the boundary 6(1) is given by 

■w»- - > " > ;V <26) 

^JJ. Mean inter-regeneration time and distribution of the content jump 

In this section we continue the analysis of the previous section and we are going to 

study the process only at the moment of the discontinuities, that is when the environment 

jumps from state 2 to state 1 and the buffer content just before the jump epoch is greater 

than the buffer level 6(1). 

These special epochs are regeneration points for the process Z(t), and we denote by 

t the general interval of time that lags betweens two such discontinuities of Z(t). We are 

going to exploit the regenerative structure of Z(t) at its discontinuity points to determine 

the expectation of r and the distribution of Z(j— ). 
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This quantities may have relevance when studying an optimization problem. For 
example, if we consider the example mentioned before, the discontinuity may model the 
start-up content of a second server to which there can be associated start-up costs taken 
into account by a given cost function. The value of 1/E[r] gives the average rate at which 
such costs incur. 

Assume that Z(O-) > 6(1) and Z(0) = 6(1) and define r = inf t>0 {AZ(t) < 0}, the 
picture below helps in understanding the given definition. 



6(2) 




Let f(z) be any twice differentiable function, by Ito's Lemma we have that 

,-t rt 



f(Z(t)) = f(Z(0)) 



a / f'(Z(s))dW(s) 



1 



<j 2 f"{Z{s))+nf{Z{s))ds 



o 



+ /'(0) L{t) - /'(&(!)) U,{t) - f{b{2)) U 2 (t) 



(27) 



where Ui(t), i G E is the amount of content lost at the barrier b(i) during the inter- 
val [0,t). Evaluating (27) at t = r under the condition that Z(0) = 6(1) and taking 
expectation we get that 



E 6(1) [/(Z(r)) -/(&(!))] = E b(1) 



j^\o 2 f"{Z{s))+»f{Z{s))ds 

+/'(0)E 6(1) [L(r)]-/'(6(l))E 6(1) [^ 1 (r)] 
-f(b(2))E b{1) [U 2 (r)] 

Having that r is a regeneration point for Z we have by renewal theory that 



(28) 



n(z) -p{z* < z} = 7 ] E b{1) 



/ l{z(s)<z} ds 
Jo 



where we defined 77 = l/£' b ( 1 )[r], the rate of discontinuities of Z. 
In general we have that 

,6(2) r ,r 

E[f(Z*)}= f(z)U(dz) = r,E H1) / f(Z(s))ds 
Jo Uo 



(29) 



so that multiplying equation (28) by 77, applying (29) and using repeatedly integration 
by parts similarly to the proof of Theorem 1, we get the following system of differential 
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equations by collecting all the integrand terms with factor f(z), 



l -a 2 tt"(z) - ^'(z) = V H'(z) l {b(1)<z < b(2)} (30) 

where H(z) = P{6(1) < Z(r-) < z}, with 6(1) < z < 6(2), and n(z) = U'(z). 

Having that U(z) = H 2 (z) + Hi(z A 6(1)), by comparing the second equation in (30) 
with the derivative of equation (23) we get that 

V H'(z)=q 21 U' 2 (z) . (31) 

Integrating last equation over the interval [6(1), 6(2)] with boundary conditions H(b(l)) = 
and H(b(2)) = 1 it follows that rj = q 2 i(ir 2 — n 2 (6(l))^. Substituting its value in (31) 
and integrating we get that, for 26/2, 

= n 2 (z)-n 2 (6(l)) = _ A(z _ b(1)) sinh(9 2 (6(2)-z)) 
{) 7r 2 -n 2 (6(l)) sinh(e 2 (6(2)-6(l))) ' 

where we used the expression of n 2 (z) given in (25). 

4-2. Two-state modulation with at least one with no diffusion component 

In this section we consider the case when N — 2 and only in one of the two states 
the diffusion coefficient is positive. 

The way to proceed to compute the stationary distribution is similar to the one used 
in the previous section, but we decided to include these examples because for these cases 
the matrices T^ 1 are rectangular and do not reduce to the identity matrix. In [4] and [5] 
the authors looked at how compute the stationary distribution for the case of two-side 
reflection with two non-modulated barriers, but no explicit examples where given and as 
far as we know they are not treated elsewhere. 

The system (10) can be solved by considering the two intervals 1\ = [0, 6(1)] and 
I 2 = [6(1), 6(2)]. We consider the two cases when the state with no diffusion component 
is the first one, i.e. o\ = and then when the state with no diffusion is the second one, 
i.e. (T2 = 0. 

In the second case in order to have positive probability for the process to enter the 
interval (6(0), 6(1)] we assume that [i 2 > 0. To simplify the analysis and reduce the cases 
we assume that the asymptotic drift k < and for the first case that (n < 0. 

Case ;«i < 0, u\ = 0; 01 > 0; n < 0: For z £ I\ we have that 

- mi ni (z) - qi2 it (z) + q 21 n 2 (z) = 

\o\ lT 2 '(z) - ^2 n 2 (z) + q 12 IT(z) - q 21 U 2 {z) = 
and for G I 2 we have 

\<*\ n 2 ( z ) - M2 n 2 (z) - q 21 n 2 (z) = -q 12 m . 

In addition we have the boundary conditions IT(0) = and IT(6(i)) = Hi — 
<Z3-m/((Zi2 +921), i — 1,2, the continuity condition II 2 (6(1)— ) = n 2 (6(l)+) and 
the differentiability condition U' 2 (b(l)-) = n' 2 (6(l)+). 

"l2 



Having \i\ < it follows that E + = {2} and E — {1,2}. Hence defining 7^ 
Pi{sup t>0 {X(t)} = 0} we have that = (7^, 1) T and T± = l~ . 

Solving the system (16) that is explicitly rewritten as 

<Zi 2 (l-7i + )-Mi Af 7 i + =0 
-?2i(l-7i + )-M2A+ + i ( r 2 2 (A+) 2 =0 
and selecting the solution with A| < we have that 



,+ _ M2 _ 912 _ 
1 ~~ ~2 ,, ,, „2 



Mi Mi 02 



1 s + 1 M2 012 2 2 °2 1 512 , 
01; 7i =-- 9-— — 01 



2 Mi 021 4q 2 iMi 2 Mi 021 



where <5i = y ,u 2 ^ 2 + 2,ui/i 2 gi 2 ^ + 4^ 2 g 2i ^ + g 2 2 ^-. Solving the system (16) 
that is explicitly rewritten as 



012 +MiA 12 

021 + M2A2J 

we get the following solutions 

012 

Mi 







io- 2 A 12 A 21 - 50 2 (A 21 ) 2 — 



M2 012 _J_ c 
A 21 = — + 1 , 2 



erj 2/xi /ii (J 2 



where h = \J 2g 2 i Mi °2 + [mi M2 + 5i2 ^o"!] ■ 

Over the interval I\ , the stationary distribution is given by 

(nM) =PlTteZXtct + PlT ' e ~ ZAT&r ■ 

Having n 2 (0) = 0we can solve for = — cj" 2 and with IIi(b(l)) = 7Ti we have 



n 2 (*) 



Pi 



n 2 (6(i)) 



(32) 



where f + = ( ° A+ = A+ I and C = [f+e b « A f - IYe"* 1 )^ 

Over the interval I 2 let A^ be the negative solutions of the equations \a\ z 2 =f^ 2 z — 
o 2 i = 0. Define 6 2 = (AJ+A2 )/2 = -( ^//x 2 + 2g 21 a- 2 )/<7 2 and A 2 = (Aj-A^)/2 - 
/U2/021 we can write the solution for z e I2 with boundary condition n 2 (6(2)) = 7r 2 
as 

In the equation above we already used the continuity condition for z = 6(1) by 
imposing the boundary value at 6(1) to be equal to n 2 (6(l)). 
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The solution is completely identified by solving for the value of 112(6(1)) that assures 
the differentiability of n 2 (z) at z = 6(1). 

We have that for z € I\ 



ti'{z) = P 1 \ttA+e zA i +IY Afe- zA ii C^P^ 1 
and for z E 7 2 



n 2 (6(i)) 



■5W - -A 2 <, - n 2(2 » + ^■''e. ^ff.-", ^ - n,«U)) . 

To simplify the resulting expression define 

k i = e 2 T P 1 [f+ A+e b(1)A ^ + 17 Afe- 6 * 1 ^] C^Pf 1 3, » = 1,2, 
where e, is the i-th column vector of the canonical base in M 2 , and 

fc 3 = e 2 coth(e 2 (6(2)-6(i))) , 

then the value of n 2 (z) at the boundary 6(1) is given by 

7T 2 (& 3 - A 2 ) - irikx 



n 2 (6(i)) 



h + h 



(34) 



Notice that in this case the function ITi(z) is not continuous at z = where it gives 
the probability to find the system empty when the environment is found in state 
1, i.e. 



Y{Z* = 0, J* = 1) = IIi(O) = e^P x 



c ^ n 



^Si)) 



Case a 1 > 0; ^ 2 > 0, er 2 = 0; n < 0: For the interval I\ we have that 

\o\ u'{(z) - mi n;(z) - gi2 ni(«) + 921 n 2 (z) = o 
-^ 2 n 2 (z) + gi2 n 1 (z)- (?21 n 2 (z)= o 

and for the interval 7 2 we have 

-,U 2 II 2 (z) - 92I II 2 (z) = -q 12 TTl . 

In addition we have the boundary conditions 1^(0) = and Hi(b(i)) = m = 
93-1,1/(912 + 921), i = 1, 2 and the continuity condition II 2 (6(1)— ) = n 2 (6(l)+). 

Having /z 2 > it follows that E + = {1,2} and E~ = {1}. In a similar way as in the 
previous case we have = 1+ and T± = (1, 7 2 ") T where 7^ = P 2 {inf t>0 {X(t)} = 
0}. Solving the system (16) that is explicitly rewritten as 



912 - Mi-^2 - \o\ A+ A+ - \a\ (A+) 



921 - M2A21 



= 
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we get the following solutions 



A+ - fl1 921 -+- 1 A • \+ - 921 

A 12 — 2 _ n 1 2 0l > A 21 — 

(T 2 2^ 2 [i 2 (j{ M2 



where <5i = \j 2q 12 \i\ a\ + [/x 2 Mi + 921 j^i] 2 . 
Solving the system (16) that is explicitly rewritten as 

-92i(l-7 2 ~) + MiAr + ^i( A r) 2 =0 

921 (1 -I2) +M2A 1 "7 2 " =0 

and selecting the solution with A^ < we have that 

x _ -Mi , 921 1 , - 1 Mi 921 92W 1 921, 

A l = T ^ 2" 2 ' T2 = ~o 1 2 + T~2 "2 

O-f fl 2 M2^f 2^2 912 4(712^2 2/4^12 



where <5 2 = yM2Mi + 2^2^1921^ +4^9i 2 ^- + 9I1 x- 
Over the interval 7i, the stationary distribution is given by 

^ n 2 (z) J -^^i 6 c i +jPll i e c i • 

Having 11(0) — we can solve for = — and knowing that III (6(1)) = m 
we get 

(n:g)=t p - <e "" Ar - e " Af > r ' <35 » 

where = I and ci = Pi (e"^ 1 ) A r - e b W A ^ )r J" . 
Over the interval I 2 the solution is simply 

n 2 (z) = 7r 2 - e -^ (z - fc(1)) (7r 2 - n 2 (6(l))) , (36) 

where n 2 (6(l)) is known by continuity from equation (35) and we used the fact 
that 7T 2 = 7ri 912/921- 

Notice that in this case the function H 2 (z) is not continuous at z = 6(2) where it 
gives the probability to find the system saturated when the environment is found 
in state 2, i.e. 

P(Z* = 6(2), J* = 2) = 7r 2 - n 2 (6(2-)) = e- q -^ {z - b{1)) (ir 2 - II 2 (6(1))) . 



5. Synchronized Barrier Strategies for Dividend Payout 

In this section we look at an application of the model presented in Section 2 to the 
problem of computing the expected dividend payouts of a company. We assume that 
the company profit fluctuates according to the Markov modulated Brownian motion 
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(X(t), J{t)) where X(t) is as in (3) and J(t) is the environment process as introduced in 
Section 2. 

The model is a direct extension of the one studied in [6] where for the profit process 
it was chosen a free Brownian motion, the introduction of the external environment J(t) 
from a modeling perspective helps in adapting the process to the case where the profit 
process may depend on external environmental situation, such as for example seasonal- 
dependent activities. 

We assume that the company select for each state of the environment J(t) a barrier 
level b(J(i)) and decides to pay dividends as soon as its surplus process reaches that 
level. In this way the surplus process behaves exactly as the Markov Modulated two- 
sided reflected Brownian motion Z(t) up to the moment of the first hitting time of the 
lower barrier, i.e. r = inf{i > : Z(t) = 0}, that corresponds to the ruin time for the 
company. 

The non-discounted total dividends paid up to time r is directly given by U(t), that 
is the upper regulator process computed at the ruin epoch. In general, given the discount 
rate 5 > 0, this amount is given by 



U is a random quantity that depends on the path realization of the process (X, J), as 
well as the selected barriers, {b(j)}j£E, and the start-up condition of the company (z,j) = 
(Z(0), J(0)). The aim of the company is to compute to then optimize the expected value 
of total discount dividend paid over its time-horizon. We denote this quantity by V(z,j), 
when then starts at time with initial capital z and in an environment state j. The 
formal definition is given by 



In the following we show how to heuristically determine the system of differential 
equation that admits as solution the desired quantity V(z,j), in the next sections we 
show how to get the same system of equations in a more rigorous way. 

Assuming that the starting position z is far from the barrier level b(J) we can assume 
that in a relatively short interval At the surplus process does not reach the reflecting 
barrier such that the following relation holds 




(37) 



V(z,j)=E[U\Z(0) 



z,J(0)=j]=E (Ztj) [U] . 



(38) 




A/. 



E {xJ) V{Z(At),J(At)) 



In addition we have that 




A/. 




fcS-E 



so that using the first order Taylor expansion of e s 



we finally get 



E (xJ) [V(Z(At), J(At))\ = V(z,j) + S V(z,j) At - ^ q jk (z - b,)+ At + o(At) 

keE 
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On the other side using the Ito formula it follows 

V(Z(At),J(At)) = V(Z(0),J(0))+n(J(0))V'(Z(0),J(0))At 
+ ^ 2 (J(0))V"{Z(0),J(0))At 

+ /Z urn v ( z (°) A K*)> fc ) A< + °( A *) 

keE 

and equating the two expressions above we get the following differential equations 
^a 2 (j)V''(z,j) + f,(j)V'(z,j)-5V(z,j) + J2q jk [V(zAb(k),k) + (z-b(k))+}=0. (39) 

k 

5.1. Regularity of "V (z , j) 

The derivation of equation (39) has been done by implicitly assuming regularity con- 
dition of the unknown function V(z,j), i.e. that it has second derivative on the interval 
(0, b(j)) with the exception of at most isolated points. In this section we prove that 
indeed V{z,j) does admits second derivative and again that it satisfies equation (39). 
In the following we assume that a(j) > as, according to (39), in the case it was zero 
we would need only the first derivative of V(z,j). The treatment below can be easily 
adapted to handle this case. 

Define T b = mm jeE {M t > {(Z(t), J(t)) = {b{j),j)}} and T = inf{i > 0,X(t) = 0} 
we have that 

V(z,j) = E {z ^le- ST » (V(b(J(T b )),J(T b )) - AX(T b )) l{T b < T }] 

with AX(t) =X(t)-X(t-). 

Define the stopping time r(h) = T z _ h AT z+h ATj with Tj = inf{t > 0, J{t-) ^ J(t)} 
and T z ±h = inf{i > 0, X(t) = z ± h}, we have that 

V(z,j) = E {Zij) [e- s ^(V(X(r(h)),J(r(h)))+AX(r(h)))} (40) 

= g + (h) V(Z + h,j) + g_(h) V(Z - h,j) + Vj{h, k) 2£ gj (h) 

keE ^ 

where 9J (h) = P (zj ){r(^) - Tj}, g±(h) = E {zJ) [e- s ^; r(h) = T z±h ], Qj = J2k^j Qjk, 
and 

Vj(h,k) = E {zj) {e- ST ->(V(X(Tj),k) + AX(Tj))\T(h) = Tj,J(Tj) = k] . 
To compute g,j{h) we have that 

gj(h) = P (ZJ) { sup \X(t) -z\<h}= ¥ {T lhl > Tj} = 1 - E [e-^ T ^} , 

0<t<T., 

where Tj/j = inf t >o{|y(i)| = h} is the hitting time of the set (— h, h) c of the process Y{t) 
that is a Brownian motion with drift fi(j) and diffusion coefficient cr(j). It is known, see 
[2], that 

E [e-^i] = cosh sech ( W^±Z^j = l - A /i 2 + o(h*) , (41) 

17 



c± cosh 



and from this it follows that gj(h) = 1 - E [e~* t h] = ^h 2 + o(h 2 ). 
To compute g±(h) we have that 

9±{h) = E {Zij) [E[e- ST ^;r(h)=T z±h \Tj}} 

POO 

= / E[e- 4T < h >; r(h) = T z±h \Tj = t] q 3 e^ 4 dt 
Jt=o 

/•OO 

= / E[e- ST ^E[l{T(h)=T z±h }\r(h),Tj = t]\Tj = t}q j e-^ t dt 
Jt=o 

and since the exit location X(T z+ h A T z _h) is independent of the exit time T z+ h A T z _h, 
see [15], assuming \i > 0, we have that 

E[l{r(/i) = T z± 4|r(/i), T 7 = t] = c± l{r(/i) < t} 

where c± = § ± 5(exp(§£ h) - l)/(cxp(§£ h) + 1), therefore 

rOO 

g±(h) = c± / E[e-* T Wl{r(/i)<t}|Tj = t]? J -e- wt dt 
Jt=o 

= c± / / e-^PlTifc, £ dij^e^** = c ± E [e- (5+ ^ )T i"i] 

"('t°5)-*(¥*/W) 

Last term in (40) is positive, it follows that 

V(z,j) < g+(h)V(z + h,j)+g-(h)V(z-h,j) 
< G+(h)V(z + h,j) + G-{h)V(z-h,j) 

where G±(h) = g±(h)/(g+(h) + g-(h)). Since G+(h) + G-(h) = 1 and G-(h) < \ < 
G + (h) the follow inequality holds 

V(z,j) - V(z + h,j) < V(z - h,j) - V(z,j) 

that implies that V(z,j) is continuous and concave in (0,b(j)), see Courant [3, Ch. IV. 2 
page. 326]. 

Rearranging terms in (40) and using that 1 — (g+(h) + g~(h)) = o(h) we have 
g_(h)[V(z,j) - V(z - h,j)} = g+(h)[V(z + h,j) - V(z,j)} + o(h) , 
and dividing it by h, letting h -> and having g±(h) — > 1/2 as h — > we finally get 

V'(z-,j) = V'(z+,j), 
i.e. the function V(z,j) is differentiable in (0,bj). Define 

lim Vi* + Kj)-2V%j) + Vi*-hJ) = ^ 
h^o h 2 
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with 

r2g+(h)-lV(z + h,j)-V(z,j) 



ip{z,j) = lim 

h—>0 



h h 

+ h h + h 2 V(Z > J) 

h 2 E fc Qjk [V(z A b kl k) + (z - b k )+] + o(h 2 ) - 
+ h 2 

where we used the fact that Vj(h, k) — > V(z A b k , k) + (z — b k ) + as h — > by bounded 
convergence. Having that (2g±(h) - l)/h ->■ ±n/a 2 and (2g_(h) + 2g_(h) - 2)/h 2 -> 
2 (5 + q 3 ) I a 2 as h -> we get that 

V>(*,j) = (z, j) - 2 (J + gj) j) + ]T gjfc [U(z A b k , k) + (z- b k )+] . (42) 

Using Schwarz's Theorem, [16], we get that V(z,j) is twice differentiate and V"(z,j) = 
ip(z,j). In addition (42) coincides with (39). 

5.2. Boundary conditions 

To determine the value of V(z,j) is necessary to add to the differential equations (39) 
two boundary conditions for each j 6 E, at the barriers and b(j). 

It is obvious that V(0, j) = 0, for any j e E and in addition we have 

V'(b(j),j) = l. 

This equation that can be found for the non-modulated case in [6] has the following 
explanation. 

Assume Az is small, starting at b(j) — Az we will touch the barrier b(j) at time T Az 
so that e~ STA * = l + o(T Az ). Defining t = TjATq, the shortest time between a change of 
state for J or the ruin epoch, we have that the process Z(t) will have the same dynamic of 
a single reflected (/x(j), Brownian motion at the barrier b(j) for < t < r that we 

denote by Y(t). The upper regulator process is known to have the following expression 

U z {t) = z + sup{0 < s < t : Y(t) V b(j) - z} 

as < t < t when Z(0) = z, notice that by definition Y(0) = 0. Since f* e - 5s U(ds) = 
U(t) +o(t), it follows that 

V(6(j),j) - V(b(j) - Az,j) = Az + E[l{M(t) < Az}(M(t) - Az)} + o(E[T Az ]), 

where M(t) = sup{0 < s < t : Y(s)} and T Az = inf{t > 0, M(t) > Az}. Having that 

-AzP{T Az < t) < E[l{M(t) < As}(M(t) - Az)] < 

we can get the following bounds 

i - nr Az > t) + °W Ta ^ < mAj)-mj)-A^1 < 1 + o(E[T Az }) ^ 

2 Az Az Az ' 

and taking the limit for Az — > and using the fact that P(Ta z > r) — > and 
o(E[T Az ]) j Az -> we obtain the result. 
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5.3. Expected dividend payout - Case \E\ = 2 

Assuming q 12 = ?2i = A, we have that for z G I\ the system of differential equations 
(39) reduces to 



±a 2 (l)V"(z, 1) + fi(l) V{(z, 1) - (A + 5)V(z, 1) + A V(z, 2) = 



\a 2 {2)V"{z, 2) + /x(2) ^'(z, 2) + A V(z, 1) - (A + 6)y(z, 2) = 



and for z e I 2 it reduces to 
1 



^ 2 (2)F"(z, 2) + M (2) V (z, 2) - (A + <5)F(z, 2) = -A(z - 6(1) + K(6(l), 1)) , 



with boundary conditions 



and regularity conditions 



1/(0,1) = V(0,2) = 
V'(b(l),l) = y'(6(2),2) = l 



V(b(l+), 2) 
V'(6(l+),2) 



y(6(i+),2) 

V'(6(l+),2). 



(43) 
(44) 



(45) 
(46) 



Let F(z) = (Xi , X 2 ) e J z (X 1 1 , —X 2 1 ) T be the matrix solution of the system of 
differential equation 



-AlF"(z) + A,F'(z) 



-(A + S) A 
A -(A + 5) 



for z € I\, with the null condition at the barrier, i.e. F(0) 
z\ < z 2 < z 3 < z 4 the four solution of the algebraic equation 



F(z) = O 

= O and defining by 



z 4 + z 3 



+ z- 
+ z 



2/i(l) , 2 M (2) 



^(1) v*(2) 
2/x(l) 2 M (2) 



4( M (1)+ M (2)) 



2 2 

+ 



(S-X)- 



a 2 (l)^(2) 



^(1) <7*(2) 

4<5(2A - 5) 



the matrix J = diag[zi, z 2 , z 3 , Z4] and the matrices Xi and X 2 are given by 



I 

V 
/ 

V 



2A/cr 2 (2) 



(21-22) (21-2:3) (Z1—Z4.) 

2(A-<5)/<t 2 (1)-z 1 (2 M (1)/<t 2 (1)+zi) 
(21 -22) (21 -23) (21 -24) 

2A/cr 2 (2) 

(23 -21) (23 -22) (23 -24) 

2(A-j)/<T 2 (l)-2 3 (2 M (l)/a 2 (l)+2 3 ) 

(23 -21) (23 -22) (23 -24) 
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2A/o- 2 (2) 



(22 -21) (22 -23) (22 -24) 

2(A-,5)/(7 2 (l)-2 2 (2 ; x(l)/<T 2 (l) + 22) 

(22 -21) (22 -23) (22 -24) 
2A/o- 2 (2) 



(24 -21) (24 -22) (24 -23) 

2(A-,5)/(7 2 (1)-24(2m(1)/<7 2 (1)+24) 

(24 -21) (24 -22) (24 -23) 



The above expressions are a rearrangement of formulas obtained by using the software 
Mathematical©. We have that 

V(z) = F(z)(k u k 2 ) T 0<z<6(l) 
V(z,2) =k 3 f(z)+g(z)-j^V(b(l),l) 6(l)<z<6(2) 

where 

and f(z) being any solution of the differential system 

^ 2 (2)/"(z) + M (2) f'(z, 2) - (A + S)f(z) = , 

with 6(1) < z < 6(2) and boundary condition /' (6(2)) = S/(X + 5). We choose the special 
solution that has 

that can be written as 

f(z) = - ^/^ cxp (6(6(2) - z)) cosh (A(6(2) - z)) , 



with0=^ and A= ^02 + ^1. 

To solve for the constants, fci, k 2 and k 3 we need to solve the following system of 
equations 

V'(b(l)-) = (1,T/'(6(1)+,2)) T 
V(b(l)-,2) = V(b(l)+,2) 

that is equal to 

F'(6(l))(fc 1 ,fc 2 ) T = (l,fc 3 .r(6(l))+.g'(6(l))) T (47) 
ejF(b(l))(h, k 2 ) T = fc 3 /(6(l))+ 3 (6(l))- x l-^ e 7F(6(l))(fc 1 ,fc 2 ) T (48) 

Using (47) we get 

(fc 1; fc 2 ) T = F'(6(l))- 1 (l,fc 3 /'(6(l)) + A/(A + ( 5)) T 
and using the (48) we get 

AM 2 ) - (1 , A + S) F(6(l)) F'(6(l))-! (A + 6 , A) T 



h = - 



(A + 6^ /(6(1)) - (A + 5) /'(&(!)) (1 , A + 5) F(6(l)) F'(6(l))-i e 2 . 
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